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Abstract 

We report several major theoretical steps towards realizing stable long-distance multichannel 
soliton transmission in Kerr nonlinear waveguide loops. We find that transmission destabilization 
in a single waveguide is caused by resonant formation of radiative sidebands and investigate the 
possibility to increase transmission stability by optimization with respect to the Kerr nonlinearity 
coefficient 7 . Moreover, we develop a general method for transmission stabilization, based on 
frequency dependent linear gain-loss in Kerr nonlinear waveguide couplers, and implement it in two- 
channel and three-channel transmission. We show that the introduction of frequency dependent loss 
leads to significant enhancement of transmission stability even for non-optimal 7 values via decay 
of radiative sidebands, which takes place as a dynamic phase transition. For waveguide couplers 
with frequency dependent linear gain-loss, we observe stable oscillations of soliton amplitudes due 
to decay and regeneration of the radiative sidebands. 

PACS numbers: 42.65.Tg, 42.81.Dp, 42.81.Qb 
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I. INTRODUCTION 


The rates of transmission of information in broadband optical waveguide systems can be 
significantly increased by transmitting many pulse sequences through the same waveguide 
13]. This is achieved by the wavelength-division-multiplexing (WDM) method, where each 
pulse sequence is characterized by the central frequency of its pulses, and is therefore called a 
frequency channel. Applications of these WDM or multichannel systems include fiber optics 
communication lines (lj-3|, data transfer between computer processors through silicon waveg¬ 
uides jd, 5], and multiwavelength lasers jfi, 7]. Since pulses from different frequency channels 
propagate with different group velocities, interchannel pulse collisions are very frequent, and 
can therefore lead to severe transmission degradation j]J]. Soliton-based transmission is con¬ 
sidered to be advantageous compared with other transmission formats, due to the stability 
and shape-preserving properties of the solitons, and as a result, has been the focus of many 
studies 1-3]. These studies have shown that effects of Kerr nonlinearity on interchannel 
collisions, such as cross-phase modulation and four-wave-mixing, are among the main im¬ 
pairments in soliton-based WDM fiber optics transmission. Furthermore, various methods 


for mitigation o 
been developed 


Kerr-induced effects, such as filtering and dispersion-management, have 


2 ), 3]. However, the problem of achieving stable long-distance propagation 
of optical solitons in multichannel Kerr nonlinear waveguide loops remains unresolved. The 
challenge in this case stems from two factors. First, any radiation emitted by the solitons 
stays in the waveguide loop, and therefore, the radiation accumulates. Second, the radiation 
emitted by solitons from a given channel at frequencies of the solitons in the other channels 
undergoes unstable growth and develops into radiative sidebands. Due to radiation accumu¬ 
lation and to the fact that the sidebands form at the frequencies of the propagating solitons 
it is very difficult to suppress the instability. I 11 the current paper, we report several major 
steps towards a solution of this important problem. 


In Refs. 


- 13j, we studied soliton propagation in Kerr nonlinear waveguide loops in the 
presence of dissipative perturbations due to delayed Raman response and nonlinear gain-loss. 
We showed that transmission stabilization can be realized at short-to-intermediate distances, 
but that at large distances, the transmission becomes unstable, and the soliton sequences 


are destroyed. Additionally, in Ref. 


y, 


we noted that destabilization is caused by reso¬ 


nant formation of radiative sidebands due to cross-phase modulation. However, the central 
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problems of quantifying the dependence of transmission stability on physical parameter val¬ 
ues and of developing general methods for transmission stabilization against Kerr-induced 
effects were not addressed. In the current paper we take on these problems for two-channel 
and three-channel transmission by performing extensive simulations with a system of cou¬ 
pled nonlinear Schrodinger (NLS) equations. We first study transmission in a single lossless 
waveguide and investigate the possibility to increase transmission stability by optimization 
with respect to the value of the Kerr nonlinearity coefficient. We then demonstrate that sig¬ 
nificant enhancement of transmission stability can be achieved in waveguide couplers with 
frequency dependent linear loss and gain and analyze the stabilizing mechanisms. This 
stabilization is realized without dispersion-management or filtering. 


II. THE COUPLED-NLS PROPAGATION MODEL 


We consider propagation of N sequences of optical pulses in an optical waveguide in the 
presence of second-order dispersion, Kerr nonlinearity, and frequency dependent linear gain- 
loss. We assume a WDM setup, where the pulses in each sequence propagate with the same 
group velocity and frequency, but where the group velocity and frequency are different for 
pulses from different sequences. The propagation is then described by the following system 
of N coupled-NLS equations QQ: 


id z ^j+d^ j +^ j \ 2 ^ j + 2'y^2\'ijj k \ 2 ^ j = iJ 7 \g j (uj)'ip j )/2, (1) 

Mi 

where i/jj is the envelope of the electric field of the jth sequence, 1 < j < N, z is propagation 
distance, t is time, u> is frequency, 7 is the Kerr nonlinearity coefficient, and the sum over 
k extends from 1 to N [l4]. In Eq. (pQ), gj(u) is the linear gain-loss experienced by the 
jth sequence, is the Fourier transform of 7/7 with respect to time, and J 7-1 is the inverse 
Fourier transform. The second term on the left-hand side of Eq. © is due to second-order 
dispersion, the third term describes self-phase modulation and intrasequence cross-phase 
modulation, while the fourth term describes intersequence cross-phase modulation. The 
term on the right-hand side of Eq. (CQ) is due to linear gain-loss. The optical pulses in the 
jth sequence are fundamental solitons of the unperturbed NLS equation id z i\)j + + 

r y\'ipj\ 2, ipj = 0. The envelopes of these solitons are given by ) = rjj exp(ixj)sech(xj), 

where xj = (t - yj - 2 faz), Xj = + fa(t ~ Vj) + ( 1 % f/2 - ^ 2 ) z, and 77 , fa, y h 
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and otj are the soliton amplitude, frequency, position, and phase. 

Notice that Eq. ([Tj) describes both propagation in a single waveguide and propagation in a 
waveguide coupler, consisting of N close waveguides 151 ]. In waveguide coupler transmission, 
each waveguide is characterized by its linear gain-loss function gj(u). The form of gj( oj) is 
chosen such that radiation emission effects are mitigated, while the soliton patterns remain 
intact. In particular, we choose the form 


9j(w) = ~9l + I (g„ + 9l) [tanh {p [w - ft(0) + W/ 2]} 

— tanh {p [u — 0j(O) — W/2]}] , (2) 


where 1 < j < N, and /3j(0) is the initial frequency of the jth sequence solitons. The 
constants gL, g e q, p, and W satisfy g^ > 0, g eq > 0, p 1, and A/3 > W > 1, where A/3 
is the intersequence frequency difference. We note that the condition A/3 > 1 is typical 


for soliton-based WDM transmission experiments 


16 


20]. Figure |T] shows typical linear 


gain-loss functions gi(u>) and g 2 ((n) for a two-channel waveguide coupler with g L = 0.5, 
g eq = 3.9 x 10 -4 , /3i(0) = —5, /3 2 (0) = 5, W = 5 and p = 10 (these parameters are used in 
the numerical simulations, whose results are shown in Fig. |8]). In the limit as p /$> 1, gj{uf) 
can be approximated by a step function, which is equal to g eq inside a frequency interval of 
width W centered about /3j( 0), and to —g^ elsewhere: 


9j{u) - 


g eq if ^(O) -W/2<co< 13,(0) + W/2, 
gL elsewhere. 


(3) 


The approximate expression (]3]) helps clarifying the advantages of using the linear gain-loss 
function fl2]) for transmission stabilization. Indeed, the relatively strong linear loss gL leads 
to efficient suppression of radiative sideband generation outside of the frequency interval 
(/3j(0) — W/2, m + W/2]. Furthermore, the relatively weak linear gain g eq in the frequency 
interval (/3j(0) — W/2,/3j(0) + W/2] compensates for the strong loss outside of this interval 
and in this manner enables soliton propagation without amplitude decay. In practice, we first 
determine the values of W, and p by performing simulations with Eqs. (P) and (J2J) with 
g eq = 0, while looking for the set that yields the longest stable propagation distance. Once 
gL, W, and p are found, we determine g eq by requiring rjj(z) = rjj( 0) = const for 1 < j < N 
throughout the propagation. More specifically, we use the adiabatic perturbation theory for 
the NLS soliton (see, e.g., Ref. [3]) to derive the following equation for the rate of change 
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FIG. 1: An example for the frequency dependent linear gain-loss functions gj(uj ) defined by Eq. 
((2|) in a two-channel waveguide coupler. The solid blue and dashed red lines correspond to g\ (cj) 
and 52 (w), respectively. 


of rjj with z due to the linear gain-loss 


drjj [ / 7 tW \ 

-£=[-9L + (s«, + SL)te*h 


r h ■ 


Requiring rjj(z) = rjj( 0) = const, we obtain the following expression for g eq : 


9eq 


tanh 


/ ttW 

V(87) 1/2 ^(0)A 


-l 


1 > 9l- 


(4) 

(5) 


Since different pnlse sequences propagate with different group velocities, the solitons un¬ 
dergo a large number of intersequence collisions. Due to the finite length of the waveguide 
and the finite separation between adjacent solitons in each sequence, the collisions are not 
completely elastic. Instead, the collisions lead to emission of continuous radiation with 
peak power that is inversely proportional to the intersequence frequency difference A/3. The 
emission of continuous radiation in multiple collisions eventually leads to pulse pattern dis¬ 
tortion and to transmission destabilization. In the current paper, we analyze the dependence 
of transmission stability on physical parameter values and develop waveguide setups, which 
lead to significant enhancement of transmission stability. 
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FIG. 2: (a) The Fourier transforms of the soliton patterns at the onset of instability \ipj(u, z u )\, 
where z u = 470, for two-channel transmission in a single lossless waveguide with 7 = 2, T = 15, 
and A/3 = 12. The blue circles and red squares represent z u )\ with j = 1,2, obtained by 

numerical solution of Eq. while the magenta diamonds and green triangles correspond to the 
theoretical prediction, (b) The soliton patterns at the onset of instability | z u ) \ for two-channel 
transmission with the parameters used in (a). The solid blue and dashed red lines correspond to 
z u )\ with j = 1 , 2 , obtained by the simulations, while the black diamonds and green triangles 
correspond to the theoretical prediction. 
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III. NUMERICAL SIMULATIONS 


Introduction 


To investigate transmission stability, we numerically integrate the system (P), using the 
split-step method with periodic boundary conditions [1]. The use of periodic boundary 
conditions means that the simulations describe pulse dynamics in a closed waveguide loop. 
The initial condition is in the form of N periodic sequences of 2/1+1 solitons with amplitudes 
r/j( 0 ), frequencies //, ( 0 ), and zero phases: 


K 


? b(0) exp[i/3j(0)(t — kT)\ 


Mt,0) cosh[( 7 /2)i/^( 0 )(( - kT )\' 


( 6 ) 


where 1 < j < N, T is the time-slot width, and N = 2 or N = 3. This initial condition 
represents the typical situation in multichannel soliton-based transmission To maxi¬ 

mize the stable propagation distance, we choose /?i ( 0 ) =-/? 2 ( 0 ) for a two-channel system, and 
/?i (0)=-/5 3 (0) and /J 2 (0) = 0 for a three-channel system. This choice is based on extensive 
numerical simulations with Eq. (JT]) with the right-hand-side set equal to zero and different 
values of /3j( 0). For concreteness, we present here the results of numerical simulations with 
parameter values T = 15, 7 /^( 0 ) = 1, K = 1, and a final transmission distance Zf = 5000. 
We emphasize, however, that similar results are obtained with other values of the physical 
parameters. That is, the results reported in this section are not very sensitive to the values 
of K, rjj( 0), and T, as long as r)j( 0) is not much smaller than or much larger than 1, and as 
long as T > 10. 


B. Two-channel transmission 

We start by considering two-channel transmission in a single lossless waveguide. Simula¬ 
tions with Eq. (JT]) with N — 2 show stable propagation at short-to-intermediate distances 
and transmission destabilization at long distances. As seen in Fig. [ 2 j the instability first 
appears as fast temporal oscillations in the soliton patterns, which is caused by resonant 
formation of radiative sidebands with frequencies /3 2 (0) for j = 1 and /3i(0) for j = 2. The 
growth of the radiative sidebands with increasing z eventually leads to the destruction of the 
soliton patterns. We note that when each soliton sequence propagates through the waveg¬ 
uide on its own, no radiative sidebands develop and no instability is observed up to distances 
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FIG. 3: (a) The Fourier transforms of the soliton patterns at the final propagation distance 
Zf )\. where Zf = 20000, in the case where each soliton sequence propagates on its own 
through a lossless waveguide. The values of the physical parameters are /3i(0) = —5, /^(O) = 5, 
7 = 2, and T = 15. The symbols are the same as in Fig. [2(a). (b) The soliton patterns at the final 
propagation distance \ipj(t,Zf)\ for the single-sequence propagation setup considered in (a). The 
symbols are the same as in Fig. [2(b). 
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FIG. 4: Stable propagation distance z s vs frequency spacing A/3 for two-channel transmission in a 
single lossless waveguide with T = 15. The black circles, orange stars, blue squares, red triangles, 
and green diamonds represent the results obtained by the simulations for 7 = 0.5, 0.75,1.0,1.5, and 
2 . 0 , respectively. 

as large as z — 20000 (see Fig. [3]). Thus, the instability is caused by the Kerr-induced in¬ 
teraction in interchannel collisions, i.e., it is associated with the intersequence cross-phase 
modulation terms in Eq. (JTj) . 

An important question about the transmission concerns the dependence of transmission 
stability on the value of the Kerr nonlinearity coefficient. In particular, we would like to find 
if there is an optimal value of 7 , which leads to minimization of radiative sideband generation 
and to maximization of transmission stability. To answer this question, we define the stable 
propagation distance z s as the distance z u at which instability develops, if z u < Zf, and as 
Zf, if no instability is observed throughout the propagation. That is, z s = z Ul if z u < z f 
(instability is observed), and z s = Zf, if z u > Zf (instability is not observed). We then 
carry out simulations with Eq. ([Tj) for N = 2, 0.5 < 7 < 2, and 4 < A/3 < 40, and plot 
z s vs frequency spacing A/3. The results of the simulations are shown in Fig. [4j It is seen 
that z s increases with increasing A/3, in accordance with the decrease of intersequence cross¬ 
phase modulation effects with increasing frequency spacing [l(J . Moreover, for all frequency 
differences A/3 in the interval 4 < A/3 < 40, the z s values obtained with 7 = 1 are larger 
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than or equal to the z s values achieved with 7 7 ^ 1. Thus, 7 = 1 is the optimal value of the 
Kerr nonlinearity coefficient. Based on these results and results of simulations with other 
sets of physical parameters, we conclude that for two-channel systems, there indeed exists 
an optimal value of the Kerr nonlinearity coefficient, which minimizes radiative sideband 
generation and yields the longest stable propagation distance. 

Since the radiative sideband for the 71 h sequence forms at frequency j3k( 0) of the other se¬ 
quence, it is very difficult to suppress radiative instability in a single waveguide by frequency 
dependent gain-loss. The situation is very different in waveguide coupler transmission, since 
in this case one can employ a different gain-loss profile for each waveguide, with strong loss 
for all frequencies outside of a frequency interval centered about /3j( 0). We therefore turn 
to consider waveguide couplers with frequency dependent linear loss, and show that in this 
case, significant enhancement of transmission stability can be achieved, even for non-optimal 
7 values. For this purpose, we numerically solve Eqs. © and ([2]) with N = 2 and g eq = 0 
for different 7 values and 4 < A/3 < 15. Here we present the results obtained for 7 = 2 , 
T = 15, gjj = 0.5, p = 10, and W = A/3/2. Similar results are obtained with other choices 
of the physical parameter values. Figure [5](a) shows the stable propagation distance z s vs 
frequency spacing A/3 as obtained in the simulations for two-channel waveguide coupler 
transmission along with the value obtained for transmission in a single lossless waveguide. 
We note that for A/3 >8 , z s — Zf — 5000, i.e. , the transmission is stable throughout the 
propagation. Moreover, the z s values obtained for waveguide coupler transmission are larger 
than the values obtained for single waveguide transmission by factors ranging between 172.2 
for A/3 = 4 and 2.22 for A/3 = 13. Additionally, as seen in Fig. [5](b), the solitons retain 
their shape throughout the propagation. 

We now turn to analyze the z dependence of soliton amplitudes for propagation in the 
waveguide coupler, since this analysis provides insight into the processes involved in trans¬ 
mission stabilization. We find three remarkably different dependences of soliton amplitudes 
on z in the frequency spacing intervals 4 < A/3 < 8, 8 < A/3 < 14, and A/3 > 14. Figure[ 6 ](a) 
shows the r)j(z) curves obtained by the simulations for three representative cases, A/3 = 4, 
A/3 = 12, and A/3 = 14. For A/3 = 4 and A/3 = 14, the soliton amplitudes decrease grad¬ 
ually to their final values. In contrast, for A/3 = 12, soliton amplitudes gradually decrease 
for 0 < z < 150, but then undergo a steep decrease in the interval 150 < z < 175, followed 
by another gradual decrease for 175 < z < 5000 [see Figures [5(a) and OF)]. To explain 
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FIG. 5: (a) Stable propagation distance z s vs frequency spacing A/3 for two-channel waveguide 
coupler transmission with frequency dependent linear loss and 7 = 2, T = 15, gL = 0.5, g eq = 0, 
p = 10, and W = A/3/2. The solid black line is the result obtained by numerical solution of 
Eqs. m and (|2j) . The dashed blue line is the result obtained by the simulations for two-channel 
transmission in a single lossless waveguide with 7 = 2 and T = 15. (b) The final pulse patterns 
\ipj(t, Zf) I, where zj = 5000, in two-channel waveguide coupler transmission with A/3 = 12. The 
symbols are the same as in Fig. [2]^b). 
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the abrupt decrease of rjj(z) in the interval 150 < z < 175, we analyze the z dependence 
of radiative sideband amplitudes, defined as |^i (/?2 ( 0 ), z)\ and |' 02 (/?i( 0 ), z)\ for j — 1 and 
j — 2 , respectively. As seen in Figures [(^c) and [3 sideband amplitudes exhibit different 
behavior for 0 < z < 120, 120 < z < 200, and 200 < z < 5000, which correspond to the 
three intervals observed for soliton amplitude dynamics. More specifically, for 0 < z < 120, 
sideband amplitudes are smaller than 10^ 3 and are slowly increasing, for 120 < z < 200 , 
sideband amplitudes increase up to a maximum of 0.321 at z = 160 and then decrease to 
below 10" 4 at z = 200 , while for 200 < z < 5000, sideband amplitudes remain smaller 
than 2 x 10” 4 . Thus, the steep drop of r) 3 (z) for 150 < z < 175 is related to the growth 
and subsequent decay of the radiative sidebands in the interval 120 < z < 200. This can 
be explained by noting that as the radiative sidebands grow, energy is rapidly transferred 
from a localized soliton form to a nonlocalized form, which is accompanied by the steep 
decay of soliton amplitudes. Additionally, the fast decay of the sidebands is a result of the 
strong linear loss Ql at frequencies /^(O) for j = 1 and /?i(0) for j = 2. Note that the sharp 
drop in rjj(z) and the associated growth and disappearance of the radiative sidebands can 
be described as a dynamic phase transition, which is similar to the transition of one phase 
of matter to another. Indeed, one can consider the solitons and the radiation to be two 
different “phases”. The abrupt disappearance of the radiation due to the presence of linear 
loss can then be viewed as a transition from an unstable transmission state, in which both 
phases exist in the waveguide, to a stable state, in which only the soliton “phase” exists. 

The waveguide couplers with net linear loss have a major disadvantage due to the decay 
of soliton amplitudes. This problem can be overcome in waveguide couplers with linear 
gain-loss by introducing the net linear gain g eq at a frequency interval of width W centered 
about (3j( 0). We investigate two-channel soliton transmission in waveguide couplers with 
linear gain-loss by performing simulations with Eqs. © and (J 2 ]) with N = 2 and g eq > 
0. To enable comparison with the results of Figures [5] and [ 6 l we discuss the results of 
simulations with the same parameter values, i.e., 7 = 2, T = 15, gL = 0.5, p = 10, and 
W = A/3/2. We find that soliton amplitudes exhibit different dynamic behavior in the 
frequency spacing intervals 4 < A/3 < 8 , 8 < A/3 < 14, and A/3 > 14, which are the same 
intervals observed for waveguide couplers with net linear loss. For 4 < A/3 < 8, amplitude 
values are approximately constant until transmission destabilization, while for A/3 > 14, 
the amplitudes are approximately constant throughout the propagation. In contrast, for 
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FIG. 6: (a) The z dependence of soliton amplitudes r/j for two-channel waveguide coupler trans¬ 
mission with frequency dependent linear loss and 7 = 2, T = 15, gL = 0.5, g eq = 0, p = 10, 
and W = A/3/2. The solid red, solid black, and dashed-dotted purple curves correspond to 771(2:) 
obtained by numerical simulations with Eqs. (Hj) and (J2J) for A/3 = 4, A/3 = 12, and A/3 = 14. The 
dashed-dotted-dotted green, circle-dashed blue, and short dashed-dotted orange curves represent 
772(2) obtained by the simulations for A/3 = 4, A/3 = 12 , and A/3 = 14. (b) Magnified versions of 
the 777(2) curves for A/3 = 12 in the interval 140 < 2 < 200. (c) The 2 dependence of radiative 
sideband amplitudes l^i(/32(0), 2)| (solid black line) and |^2(/3 i( 0), z )\ (dashed blue line), obtained 
by the simulations for A/3 = 12. 
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FIG. 7: The Fourier transforms of the soliton patterns \ipj(cj,z)\ in two-channel waveguide coupler 
transmission with frequency dependent linear loss for A/3 = 12 and the same values of 7, T, gL, 
g eq , p, and W as in Fig. [6l (a) \4>j(<jj,z)\ at z = 160. (b) \ipj(uj,z)\ at z = 225. The symbols are 
the same as in Fig. Oja) 

8 < A/3 < 14, the amplitudes exhibit stable oscillations throughout the propagation. Figure 
ETa) shows the oscillatory dynamics for A/3 = 10 and g eq = 3.9 x 10 4 . As can be seen, 
the amplitudes undergo a steep decrease, followed by oscillations about the value r] s = 0.86. 
Additionally, as seen in Fig. OJb), pulse distortion at Zf = 5000 is small, although the 
solitons within each sequence experience position shifts relative to one another. To check if 
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the oscillations of soliton amplitudes are caused by radiative sideband dynamics, we analyze 
the 2 dependence of radiative sideband amplitudes l^i(/?2(0), z)\ and |'02(/3 i(O), z)\. As seen 
in Fig. |8|(c), the amplitudes of the radiative sidebands experience alternating “periods” of 
growth and decay. Furthermore, the points where the sidebands are maximal are located 
near the beginnings of the relatively short intervals, where soliton amplitudes are decreasing 
[see Fig. [81(a)]. Based on these observations, we conclude that the oscillatory dynamics of 
soliton amplitudes is caused by decay and regeneration of the radiative sidebands. This can 
be explained by noting that as the sidebands grow, energy is transferred from a localized 
soliton form to a nonlocalized form. The strong linear loss ()j J outside the central frequency 
intervals leads to relatively fast decay of the radiative sidebands, which is accompanied 
by a decrease in soliton amplitudes. Furthermore, the weak linear gain g eq at the central 
frequency intervals leads to slow growth of soliton amplitudes at the subsequent waveguide 
spans and to the observed oscillatory dynamics. 

C. Three-channel transmission 

It is important to investigate whether the results obtained in subsection IIIIBI for trans¬ 
mission stabilization in a two-channel system remain valid as more frequency channels are 
added. For this purpose, we turn to discuss the results of numerical simulations for three- 
channel transmission, starting with transmission in a single lossless waveguide. As seen in 
Fig. [9l transmission destabilization is caused by resonant formation of radiative sidebands 
in a manner similar to the two-channel case. Moreover, the largest radiative sidebands of 
the jth sequence appear at frequencies /?*,(0) of the neighboring soliton sequences. That is, 
the largest sideband of the j — 1 sequence is formed at frequency /3 2 (0), the j — 2 sidebands 
are formed at frequencies /?i(0) and /3 3 (0), and the j = 3 sideband is formed at frequency 
/3 2 (0). Similar to the two-channel case, the formation of the radiative sidebands leads to 
pulse distortion, which first appears as fast oscillations in the soliton patterns. The growth 
of the radiative sidebands with increasing propagation distance eventually leads to the de¬ 
struction of the soliton patterns. Furthermore, the distances z u , at which instability first 
appears in three-channel transmission, are significantly shorter compared with the distances 
z u observed for two-channel transmission. For example, for parameter values 7 = 2, T = 15, 
and A (3 = 12, used in Figs. [2] and [9] z u = 74 for A = 3 compared with z u = 470 for N = 2. 
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FIG. 8: (a) The z dependence of soliton amplitudes iy for two-channel waveguide coupler trans¬ 
mission with frequency dependent linear gain-loss and 7 = 2, T = 15, A f3 = 10, 57 = 0.5, 
g eq = 3.9 x 10 , p = 10, and W = 5. The solid black and dashed blue lines correspond to rjj(z) 

with j = 1,2, as obtained by numerical solution of Eqs. © and ({2j) . The green circles indicate the 
distances at which radiative sideband amplitudes attain their maxima, (b) The final pulse patterns 
Zf)\, where Zf = 5000. The symbols are the same as in Fig. [2jb). (c) The z dependence of 
radiative sideband amplitudes |V’i(/?2(0), z)\ (solid black line) and IV^C/MO), z)\ (dashed blue line), 
obtained by the simulations. 
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FIG. 9: (a) The Fourier transforms of the soliton patterns at the onset of instability z u )\. 

where z u = 74, for three-channel transmission in a single lossless waveguide with 7 = 2, T = 15, 
and A/3 = 12. The blue circles, red squares, and orange stars represent z u )\ with j = 1,2, 

and 3, obtained by numerical solution of Eq. ©, while the magenta diamonds, green triangles 
and black crosses correspond to the theoretical prediction, (b) The soliton patterns at the onset of 
instability | z u )\ for three-channel transmission with the same parameters used in (a). The solid 
blue, dashed red, and dash-dot purple lines correspond to z u )\ with j = 1,2, and 3, obtained 
by the simulations, while the black diamonds, green triangles and orange stars correspond to the 
theoretical prediction. 
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FIG. 10: Stable propagation distance z s vs frequency spacing A (3 for three-channel transmission in 
a single lossless waveguide with T = 15. The black circles, orange stars, blue squares, red triangles, 
and green diamonds represent the results obtained by the simulations for 7 = 0.5, 0.75,1.0,1.5, and 
2 . 0 , respectively. 

Next, we discuss the dependence of transmission stability for three-channel transmission 
in a single lossless waveguide on the frequency spacing A/3 and the Kerr nonlinearity coef¬ 
ficient 7 . Figure fTTTI shows the stable propagation distances z s as functions of the frequency 
spacing A/3 for T — 15 and 7 = 0.5,0.75,1.0,1.5, 2.0. It is observed that the largest z s 
values are obtained with 7 = 0.5 for 4 < A/3 < 6 , 11 < A/3 < 13, and 25 < A/3 < 40; 
with 7 = 1.0 for 6 < Af3 < 11, 13 < Af3 < 15, and 17 < A/3 < 25; and with 7 = 1.5 for 
15 < A/3 < 17. Thus, there is no single value of 7 , which is optimal over the entire frequency 
spacing interval 4 < A (3 < 40. This behavior is sharply different from the one observed for 
two-channel transmission, where the value 7 = 1.0 is found to be optimal over the entire 
interval 4 < A [3 < 40. Furthermore, the z s values obtained for N = 3 are significantly 
smaller than the ones obtained for N = 2. For example, for 20 < A/3 < 40 the z s values for 
three-channel transmission are smaller than 1000 for all 7 values, while the corresponding 
z s values for two-channel transmission are equal to 5000 for 7 = 1.0,1.5, and 2.0. 

We now turn to analyze three-channel transmission in a waveguide coupler with frequency 
dependent linear loss. Our goal is to check whether the introduction of frequency dependent 
linear loss in a waveguide coupler leads to enhancement of transmission stability in three- 
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channel systems. For this purpose, we numerically solve Eqs. © and © with N = 3, 
g eq = 0, and gi = 0.5 for 4 < A/3 < 40. To enable comparison with the results obtained for 
two-channel transmission, we present here the results of simulations with the same physical 
parameter values as the ones used in Figs. Eland EH That is, we use 7 = 2, T = 15, p = 10, 
and W = A/3/2. Figure ITTT a) shows the stable propagation distance z s vs frequency spacing 
A/3 as obtained in the simulations together with the values obtained for transmission in a 
single lossless waveguide. We observe that z s = Zf = 5000 for all A/3 values in the interval 
4 < A/3 < 40. Additionally, as seen in Fig. [TTT b). pulse-pattern distortion is relatively 
small at the final propagation distance. Based on these observations we conclude that three- 
channel transmission through the waveguide coupler is stable throughout the propagation 
for any A/3 value in the interval 4 < A/3 < 40. Surprisingly, the z s values obtained for three- 
channel waveguide coupler transmission for 4 < A/3 < 7 are larger than the corresponding 
values obtained for two-channel waveguide coupler transmission by factors ranging between 
6.45 for A/3 = 4 and 1.25 for A/3 = 7. Furthermore, the z s values obtained for three-channel 
waveguide coupler transmission are larger than the values obtained for three-channel single 
waveguide transmission by factors ranging between 1250 for A/3 = 4 and 9.43 for A/3 = 30. 
Note that these enhancement factors are significantly larger than the enhancement factors 
for two-channel transmission, which are smaller than 172.3 for all A/3 values in the interval 
4 < A/3 < 40. 

Further insight into the mechanisms leading to transmission stabilization in waveguide 
couplers with frequency dependent linear loss is gained by analyzing the z dependence of 
soliton amplitudes. Similar to the two-channel case, we find three qualitatively different 
dependences of soliton amplitudes on 2 in the frequency spacing intervals 4 < A/3 < 8 , 
8 < A/3 < 14, and A/3 > 14. FigureEHa) shows the r)j(z) curves obtained by the simulations 
for three representative A/3 values, A/3 = 4, A/3 = 12 , and A/3 = 14. We observe that for 
A/3 = 4 and A/3 = 14, the soliton amplitudes gradually decrease to their final values. 
For A/3 = 12, the amplitudes of the solitons in the first frequency channel also decrease 
gradually throughout the propagation. However, the amplitudes of the solitons in the second 
and third frequency channels exhibit a more complicated dependence on z, which is very 
similar to the one observed for two-channel waveguide coupler transmission with A/3 = 12 . 
More specifically, soliton amplitudes in the second and third channels gradually decrease for 
0 < z < 150, but then undergo a steep decrease in the interval 150 < z < 200, followed by 
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FIG. 11: (a) Stable propagation distance z s vs frequency spacing A/3 for three-channel waveguide 
coupler transmission with frequency dependent linear loss and 7 = 2, T = 15, 57 = 0.5, g eq = 0, 
p = 10, and W = A/3/2. The solid black line is the result obtained by numerical solution of 
Eqs. © and ©. The dashed blue line is the result obtained by the simulations for three-channel 
transmission in a single lossless waveguide with 7 = 2 and T = 15. (b) The final pulse patterns 
\ipj(t, Zf)\, where Zf = 5000, in three-channel waveguide coupler transmission with A/3 = 12. The 
symbols are the same as in Fig. [9l/b). 
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another gradual decrease for 200 < z < 5000 [see Figures fl2faf and fl2f b)]. 

The behavior of r/j(z) in the interval 150 < z < 200 can be explained by analyzing the 
z dependence of radiative sideband amplitudes |^(/3jk(0), z)\, where 1 < j < 3, 1 < k < 3, 
and j 7 ^ k. Figure fTTf c) shows the z dependence of the radiative sideband amplitudes in the 
interval 125 < z < 250, while Fig. IT31 shows the Fourier transforms of the soliton patterns 
at z = 175 and z = 250. As can be seen from these figures, the sideband amplitudes 
1^2 (^ 3 ( 0 ), z)\ and |^ 3 (/ 02 (O), z)\ attain a sharp maximum at z — 177.5 with maximal values 
of 0.456 and 0.465, respectively. The increase of these sideband amplitudes is followed by a 
drop to values smaller than 10 -3 at z — 205. The formation and subsequent decay of the 
main radiative sidebands for the j = 2 and j = 3 channels in the interval 150 < z < 200 
explains the sharp drop in rj- 2 (z) and r/ 3 (z) observed in this interval. Indeed, the formation 
of the sidebands leads to energy transfer from a localized soliton form to a nonlocalized 
radiative form, which results in the steep drop of r) 2 (z) and rj 3 (z). Additionally, the strong 
linear loss ()l at frequencies /^(O) for j = 2 and /3 2 ( 0) for j = 3 leads to the relatively fast 
decay of the sidebands following their formation. We note that the evolution of rj 2 (z) and 
r) 3 (z) in the three-channel waveguide coupler is in fact quite similar to the evolution of rfr(z) 
and rj 2 (z) in the two-channel waveguide coupler. 

Figure fl2fcl also shows that the sideband amplitudes |Vh(/M0), z)\ and |' 02 (/?i(O), z)\ 
attain a maximum at z — 185 with maximal values of 0.117 and 0.111, respectively. The 
increase of these sideband amplitudes is followed by a decrease to below 10~ 3 values at 
z = 257.5. Thus, the formation and subsequent decay of the j = 1 sideband at frequency 
^ 2 ( 0 ) in an interval centered about z = 185 explains the observed drop in the value of rji(z) 
in this interval. Additionally, the sideband amplitudes |^i(/3 3 (0), z)| and |t/b(/M0), z) \ attain 
a maximum at z = 177.5, but the corresponding maximal values are smaller than 10” 3 , and 
as a result, do not significantly affect the amplitude dynamics. The relatively small values of 
l^i(£ 3 ( 0 ), 2 )| and | , 03 (/ 3 i ( 0 ), z)\ compared with the other four sideband amplitudes indicate 
that the magnitude of radiative sidebands decreases as the absolute value of the frequency 
difference |/3j(0) — ^( 0)1 increases. 

We conclude the discussion of three-channel transmission by considering propagation in 
waveguide couplers with frequency dependent linear gain and loss. As explained in section 
m in these waveguide couplers, the weak linear gain g eq in the frequency interval (/7/(0) — 
VF/2,/?j(0) + Wj 2] is expected to enable soliton propagation without amplitude decay. To 
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FIG. 12: (a) The z dependence of soliton amplitudes r/j for three-channel waveguide coupler 

transmission with frequency dependent linear loss and 7 = 2, T = 15, 57 = 0.5, g eq = 0, p = 10, 
and W = A/3/2. The solid red, dashed black, and solid purple curves correspond to 771 ( 2 ) obtained 
by numerical simulations with Eqs. (HJ and © for A/3 = 4, A/3 = 12 , and A/3 = 14. The 
dashed-dotted-dotted green, short dashed blue, and short dashed-dotted orange curves represent 
772 ( 2 :) obtained by the simulations for A/3 = 4, A/3 = 12 , and A/3 = 14. The dotted blue, dashed- 
dotted purple, and dotted green curves represent 773 ( 2 :) obtained by the simulations for A/3 = 4, 
A/3 = 12 , and A/3 = 14. (b) Magnified versions of the 77 (z) curves for A/3 = 12 in the interval 
125 < z < 250. (c) The z dependence of radiative sideband amplitudes |V’i (/32 (0), ^r) | (solid blue 
line), |t/>i(/ 33 ( 0 ), z)\ (solid purple line), |V , 2 (/3i( 0 ) , 2 :)| (dashed red line), IV^C^O), z)\ (solid black 
line), 1 7/3 (/3i(0), 2 ) | (dashed light blue line), and | (^ 2 ( 0 ), z) | (dashed green line) obtained by the 
simulations for A/3 = 12 . 
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FIG. 13: The Fourier transforms of the soliton patterns \ipj(u,z)\ in three-channel waveguide 
coupler transmission with frequency dependent linear loss for A/3 = 12 and the same values of 7 , 
T, gLt 9eqi Pi and W as in Fig. [12j (a) \'tpj(oj,z)\ at z = 175. (b) \^>j(u,z)\ at z = 250. The 
symbols are the same as in Fig. [9]^a) 

check if such stable propagation can indeed be realized, we numerically solve Eqs. (P 
and (J2J) with N = 3 and with a value of g eq , which is determined by Eq. ((HJ). To enable 
comparison with the results presented in Fig. [H] for two-channel transmission, we discuss 
the results of numerical simulations with the same set of physical parameter values. That 
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is, we use 7 = 2, T = 15, A/3 = 10, g eq = 3.9 x 10 4 , 57 , = 0.5, p = 10, and W = 5. Figure 
da) shows the z dependence of soliton amplitudes obtained by the simulations. It is seen 
that the amplitudes undergo a sharp drop, which is followed by oscillations about values 
of rjsi — 0.940, r ] S 2 = 0.768, and g s3 = 0.947, for j — 1, j — 2 and j = 3, respectively. In 
addition, as seen in Fig. [141( b). the soliton shape is retained at Zf = 5000, although the 
pulses in each sequence experience significant position shifts relative to one another. We 
note that rj S 2 is significantly smaller than r] s 1 and rj s 3 . I 11 addition, the overall oscillatory 
dynamics of soliton amplitudes is similar to the one observed in two-channel transmission, 
although the pattern of oscillations is more complex for N — 3 compared with N — 2. 

Similar to the situation in two-channel transmission, the oscillations of soliton amplitudes 
can be related to radiative sideband dynamics. We study this dynamics by analyzing the 
^-dependence of sideband amplitudes \'i/jj(/3k(0), z)\, where 1 < j < 3, 1 < k < 3, and 
j 7 ^ k. Figure |T5]( a) shows the z dependence of sideband amplitudes |'0i( / d 2 (O), z)| and 
|V’2(/di(0),z)|, Fig. fl5fb) shows the z dependence of ^(/^(O), z)| and |'03(/3 2 (O), z)|, while 
Fig. [151( c) shows the z dependence of |-0 i(/3 3 (O), z)| and |'03(/h(O), z)\. All curves in Fig. 
m are obtained by numerical solution of Eqs. (jTj) and (J2)) . We note that the values of 
|'0 i(/3 3 (O), z) I and |-0 3 (/3i(O), z)| are smaller than 0.041 throughout the propagation, and 
therefore these sidebands do not significantly affect amplitude dynamics. We therefore 
focus attention on dynamics of the four strongest sidebands |-0i(/3 2 (0), z)|, |' 02 (/?i(0), z)|, 
|'0 2 (/3 3 (O),z)|, and |'0 3 (/3 2 (0), z)|. As seen in Figure HH the radiative sideband amplitudes 
experience alternating “periods” of growth and decay, similar to the situation in two-channel 
transmission. Furthermore, the distances at which sideband amplitudes attain their maxima 
for the four strongest sidebands are located inside the relatively short intervals, where soliton 
amplitudes are decreasing. Therefore, the dynamics of the radiative sidebands can indeed be 
related to the oscillatory dynamics of soliton amplitudes. More specifically, as the sidebands 
grow, energy is transferred from a localized soliton form to a nonlocaized form, resulting 
in a decrease in soliton amplitudes. The strong linear loss 57 outside the central frequency 
intervals leads to a relatively fast decay of the sidebands, and as a result, the sidebands 
maxima are very narrow with respect to z. Additionally, the weak linear gain g eq at the 
central frequency intervals leads to the slow growth of soliton amplitude at the subsequent 
waveguide spans and to the overall oscillatory dynamics. 

Figure [141 also provides an explanation for the smaller value of rj S 2 compared to and 
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FIG. 14: (a) The z dependence of soliton amplitudes gj for three-channel waveguide coupler trans¬ 
mission with frequency dependent linear gain-loss and 7 = 2, T = 15, A/3 = 10, gL = 0.5, 
g eq = 3.9 x 10 , p = 10, and W = 5. The solid black, dashed blue, and dashed-dotted red lines 
correspond to rjj(z) with j = 1,2,3, as obtained by numerical solution of Eqs. P) and (J2J). (b) 
The final pulse patterns Zf)\, where Zf = 5000. The symbols are the same as in Fig. [9)/b). 

77 , 53 . Indeed, during the first (and largest) drop in soliton amplitudes, the solitons in the 
j = 2 sequence lose energy due to formation of radiative sidebands at both /3i(0) and / 3 3 ( 0 ). 
In contrast, the j = 1 and j = 3 solitons lose energy almost entirely due to formation 
of radiative sidebands at /3 2 (0), since the sidebands at /? 3 (0) for j = 1 and at /?i(0) for 
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FIG. 15: The z dependence of radiative sideband amplitudes for three-channel waveguide cou¬ 
pler transmission with the same physical parameter values as in Fig. [TTJ (a) |Vh(/32(0), z)\ and 
|V'2(/3i(0),z)| vs z. (b) |V> 2 (/33(0 ),z)| and 1^3(#2(0), z)\ vs z. (c) |'0i(^ 3 (O), z)\ and \$ 3 (Pi(0), z)\ vs 
z. All curves represent results obtained by simulations with Eqs. dTJ) and (j2|) . The symbols are the 
same as in Fig. He). 

j — 3 are very small. In addition, the complicated pattern of radiative sideband growth and 
decay shown in Fig. [15] is responsible for the more complex pattern of amplitude oscillations 
observed in three-channel transmission compared with two-channel transmission [compare 
Fig. Ha) with Fig. [HX a ) and Fig. IT51 with Fig. [S)[c)]. 
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IV. CONCLUSIONS 


In summary, we made several major theoretical steps towards realizing stable long- 
distance multichannel soliton transmission in Kerr nonlinear waveguide loops. We found 
that transmission destabilization in a single lossless waveguide is caused by resonant forma¬ 
tion of radiative sidebands clue to intersequence cross-phase modulation. We then showed 
that in two-channel systems, significant enhancement of the stable propagation distance, 
which holds over a wide range of interchannel frequency spacing values, is obtained by op¬ 
timization with respect to the Kerr nonlinearity coefficient 7 . In contrast, we found that 
in three-channel transmission in a single lossless waveguide, no single value of the Kerr 
nonlinearity coefficient is optimal for the entire interval of interchannel frequency spacings 
that we examined. Moreover, we developed a general method for transmission stabilization, 
based on frequency dependent linear gain-loss in Kerr nonlinear waveguide couplers, and 
implemented the method in two-channel and three-channel transmission. We showed that 
the introduction of frequency dependent loss leads to significant enhancement of transmis¬ 
sion stability even for non-optimal 7 values via decay of radiative sidebands, which can 
be described as a dynamic phase transition. For waveguide couplers with frequency de¬ 
pendent linear gain-loss, we observed stable oscillations of soliton amplitudes due to decay 
and regeneration of radiative sidebands. Transmission stabilization was achieved without 
dispersion-management or filtering. 
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